Search for: All records

Abstract Accurate pulsar astrometric estimates are essential to almost all high-precision pulsar timing experiments. Traditional pulsar timing techniques refine these estimates by including them as free parameters when fitting a model to observed pulse time-of-arrival measurements. However, reliable submilliarcsecond astrometric estimations require years of observations. Even then, power from red noise can be inadvertently absorbed into astrometric parameter fits. This effect biases the resulting estimates and reduces the sensitivity to red noise processes, including gravitational waves (GWs). In this work, we seek to mitigate these shortcomings by using pulsar astrometric estimates derived from very long baseline interferometry (VLBI) as priors for the timing fit. First, we used VLBI and timing astrometric estimates of 18 millisecond pulsars to calibrate a rotation between the reference frames used in timing and VLBI, with a precision of ∼0.7 mas. Through this frame tie, we combined timing- and VLBI-based probabilities to obtain a maximum-posterior astrometric solution. We found offsets between our results and the timing-based astrometric solutions, which, if real, would lead to the absorption of spectral power at the frequencies of interest for single-source GW searches. However, we do not find significant power absorption due to astrometric fitting at the low-frequency domain of the GW background. 

Abstract We simulate scattering delays from the interstellar medium to examine the effectiveness of three estimators in recovering these delays in pulsar timing data. Two of these estimators use the more traditional process of fitting autocorrelation functions to pulsar dynamic spectra to extract scintillation bandwidths, while the third estimator uses the newer technique of cyclic spectroscopy on baseband pulsar data to recover the interstellar medium’s impulse response function. We find that either fitting a Lorentzian or Gaussian distribution to an autocorrelation function or recovering the impulse response function from the cyclic spectrum are, on average, accurate in recovering scattering delays, although autocorrelation function estimators have a large variance, even at high signal-to-noise ratio (S/N). We find that, given sufficient S/N, cyclic spectroscopy is more accurate than both Gaussian and Lorentzian fitting for recovering scattering delays at specific epochs, suggesting that cyclic spectroscopy is a superior method for scattering estimation in high-quality data. 

Abstract PINTis a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework withinPINTto characterize the single-pulsar noise processes present in pulsar timing data sets. This framework enables parameter estimation for both uncorrelated and correlated noise processes, as well as model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated data sets as well as a real data set of PSR B1855+09. We also describe the new features implemented inPINTsince it was first described in the literature. 

Abstract Free-floating objects (FFOs) in interstellar space—rogue planets, brown dwarfs, and large asteroids that are not gravitationally bound to any star—are expected to be ubiquitous throughout the Milky Way. Recent microlensing surveys have discovered several free-floating planets that are not bound to any known stellar systems. Additionally, three interstellar objects, namely 1I/’Oumuamua, 2I/Borisov, and 3I/ATLAS, have been detected passing through our solar system on hyperbolic trajectories. In this work, we search for FFOs on hyperbolic orbits that pass near millisecond pulsars (MSPs), where their gravitational influence can induce detectable perturbations in pulse arrival times. Using the NANOGrav 15 yr narrow band dataset, which contains high-precision timing data for 68 MSPs, we conduct a search for such hyperbolic scattering events between FFOs and pulsars. Although no statistically significant events were detected, this nondetection enables us to place upper limits (ULs) on the number density of FFOs as a function of their mass within our local region of the Galaxy. For example, the UL on the number density for Jupiter-mass FFOs (∼10^−2.5–10^−3.5M_⊙) obtained from different pulsars ranges from 5.25 × 10⁶pc⁻³to 5.37 × 10⁹pc⁻³, while the UL calculated by combining results from all the pulsars is 6.03 × 10⁵pc⁻³. These results represent the first constraints on FFO population derived from pulsar timing data. 

Abstract Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 ± 0.8 mm s^–1yr^–1, which deviates from the acceleration predicted by equilibrium Milky Way (MW) models at a 4σlevel. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the MW interacting with orbiting dwarf galaxies (∼1 mm s^–1yr^–1), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in Gaia DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 60 kyr, a companion mass of 0.3M_⊙(indicative of a white dwarf or low-mass star), and a semimajor axis of 1900 au. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios. 

Abstract Pulse profile stability is a central assumption of standard pulsar timing methods. Thus, it is important for pulsar timing array experiments such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) to account for any pulse profile variability present in their data sets. We show that in the NANOGrav 15 yr data set, the integrated pulse profile of PSR J1022+1001 as seen by the Arecibo radio telescope at 430, 1380, and 2030 MHz varies considerably in its shape from observation to observation. We investigate the possibility that this is due to the “ideal feed assumption” (IFA), on which NANOGrav’s routine polarization calibration procedure relies. PSR J1022+1001 is  ∼90% polarized in one pulse profile component, and also has significant levels of circular polarization. Time-dependent deviations in the feed’s polarimetric response (PR) could cause mixing between the intensityIand the other Stokes parameters, leading to the observed variability. We calibrate the PR using a mixture of measurement equation modeling and measurement equation template matching techniques. The resulting profiles are no less variable than those calibrated using the IFA method, nor do they provide an improvement in the timing quality of this pulsar. We observe the pulse shape in 25 MHz bandwidths to vary consistently across the band, which cannot be explained by interstellar scintillation in combination with profile evolution with frequency. Instead, we favor phenomena intrinsic to the pulsar as the cause. 

Abstract We test the impact of an evolving supermassive black hole mass scaling relation (M_BH–M_bulge) on the predictions for the gravitational-wave background (GWB). The observed GWB amplitude is 2–3 times higher than predicted by astrophysically informed models, which suggests the need to revise the assumptions in those models. We compare a semi-analytic model’s ability to reproduce the observed GWB spectrum with a static versus evolving-amplitudeM_BH–M_bulgerelation. We additionally consider the influence of the choice of galaxy stellar mass function (GSMF) on the modeled GWB spectra. Our models are able to reproduce the GWB amplitude with either a large number density of massive galaxies or a positively evolvingM_BH–M_bulgeamplitude (i.e., theM_BH/M_bulgeratio was higher in the past). If we assume that theM_BH–M_bulgeamplitude does not evolve, our models require a GSMF that implies an undetected population of massive galaxies (M_⋆≥ 10¹¹M_⊙atz> 1). When theM_BH–M_bulgeamplitude is allowed to evolve, we can model the GWB spectrum with all fiducial values and anM_BH–M_bulgeamplitude that evolves asα(z) =α₀(1 +z)^1.04±0.5. 

Abstract We present the first targeted searches for continuous gravitational waves (CWs) from 114 active galactic nuclei that may host supermassive black hole binaries, using the NANOGrav 15 yr dataset. By incorporating electromagnetic priors on sky location, distance, redshift, and CW frequency, our strain and chirp-mass upper limits are typically improved by a factor of ∼2 (median 2.2) relative to all-sky limits at the same frequency. Bayesian comparisons against a model including only a Hellings–Downs-correlated background disfavors a CW signal for all targets, with a mean Bayes factor of 0.73 ± 0.32. Two targets have Bayes factors slightly above unity, but coherence tests, random-targeting experiments, and a conservative accounting of the 114-target trials factor all indicate that they are consistent with noise. We use these two candidates as worked examples to illustrate an end-to-end targeted CW search analysis and a suite of follow-up tests that future promising candidates would need to pass. We find that the electromagnetic interpretations of both candidates are ambiguous, and we update the constraints on a putative binary in 3C 66B, ruling out part of its previously allowed parameter space. Ultimately, our results demonstrate the current sensitivity of targeted pulsar timing array searches for CWs and define a road map for future multimessenger CW detections.